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Motivated by experiments exploring the physics of neutral atoms in artificial magnetic fields, 
we study the ground state of bosons interacting with long range dipolar interactions on a two-leg 
ladder. We focus on the limit where the number of particles per site is large, and the interactions 
are weak. Using two complimentary variational approaches, we find rich physics driven by the long 
range forces. Generically, long range interactions tend to destroy the Meissner phase in favor of 
modulated density wave phases. For example, nearest neighbor interactions produce an interleg 
charge density wave phase, where the total density remains uniform, but the density on each leg of 
the ladder is modulating in space, out-of-phase with one another. Next nearest neighbor interactions 
lead to a fully modulated biased ladder phase, where all the particles are on one leg of the ladder. 

This state simultaneously breaks Z -2 reflection symmetry and U{1) translation symmetry. For values 
of the flux near cf = ir, we find a switching effect for arbitrarily weak interactions, where the density 
is modulated, but the chiral current changes sign on every plaquette. Arbitrarily weak attractive 
interactions along the rungs destroy the Meissner phase completely, in favor of a modulated density 
wave phase. Remarkably, varying the rung to ladder hopping produces a cascade of first order 
transitions between modulated density wave states with different wave-vectors, which manifests 
itself as discrete jumps in the chiral current. Polarizing the dipoles along the ladder direction yields 
a region of phase space where a stable biased ladder phase occurs even at arbitrarily weak rung 
hopping. We discuss experimental consequences of our work, and draw connections between our 
work and recent experiments on cold atoms in synthetic dimensions. 


I. INTRODUCTION 

The study of interacting bosons under rotation is a 
fundamental and rich problem: for weak rotation and 
weak interactions, one finds a rich array of vortex phases, 
whereas at strong rotation and strong interactions, one 
obtains bosonic analogues of the fractional quantum Hall 
effect [I|. More recently, it has been realized that syn¬ 
thetic magnetic fields or spin-orbit coupling generated by 
coupling atoms to Raman lasers, can mimic rotation, in 
that it can lead to “kinetic frustration” by flattening the 
single particle band or introducing multiple degenerate 
minima in the band structure (2l-l9|. The presence of a 
large single-particle degeneracy implies that interaction 
effects are crucial in determining the ground state (see 
Ref. and references therein). The recent cooling 

of polar molecules [HHIl, magnetic and Ryd¬ 

berg atoms mill now offers the unique opportunity to 
explore the interplay between novel single-particle band 
structures and long range interactions in systems that 
have no traditional solid state analogues. Here we study 
the simplest example of such a system, which can be 
readily realized in current experiments: a spinless dipolar 
Bose gas trapped on a two-leg ladder in a large magnetic 
field. 

Motivated by ongoing experiments on polar molecules, 
Rydberg and magnetic atoms, the rich physics of dipolar 
gases in low dimensional systems has recently been ex¬ 
plored by several authors: in arrays of one-dimensional 
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(ID) optical lattices, at low densities, strong density- 
density interactions give rise to ordered crystalline phases 
at rational filling fractions forming a devil’s staircase 
[TqI - I^ . By taking advantage of the low-lying rota¬ 
tional states of polar molecules, and the anisotropy of 
the dipolar interaction, spin Hamiltonians such as the 
XXZ spin chain with direction dependent couplings, or 
the bi-quadratic spin-1 Haldane Hamiltonian can be real¬ 
ized in deep lattices, allowing the study of symmetry pro- 
tected topological phases on two leg ladders . In 

engineered lattice potentials with non-triyial band topol¬ 
ogy, dipolar interactions naturally giye rise to lattice ana¬ 
logues of the fractional quantum Hall effect, namely frac¬ 
tional Chern insulators (^ . In the continuum, spin- 
orbit coupled bosons with long range interactions pos¬ 
sess ground states with noyel topological defects or quasi¬ 
crystalline order [2^, [2^ . Our study complements these 
earlier works by focussing on the weak coupling limit of 
large onsite occupation, where the interplay between ki¬ 
netic frustration and long range interactions on a two 
leg ladder leads to superfluid phases with broken trans¬ 
lational and reflection symmetries. 

Experimentalists at Munich recently engineered a sys¬ 
tem of bosons on a two-leg ladder by using Raman lasers 
to create a uniform magnetic flux (f> per plaquette Q. 
They explored the non-interacting physics at a fixed flux 
(j) = TT/ 2 and found two phases as a function of the rung- 
to-leg coupling strength: a saturated chiral current or 
Meissner phase at large coupling, where equal and oppo¬ 
site currents flow on each leg of the ladder, and a modu¬ 
lated density or yortex phase at small coupling, where the 
density is modulated along the legs of the ladder. These 
phases were first discussed theoretically by Orignac and 
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Giamarchi (2^ using bosonization, although their exis¬ 
tence had already been predicted in the context of weakly 
coupled Josephson junction arrays [^ . The strongly in¬ 
teracting limit of this problem has now been thoroughly 
explored, where the interplay between single-particle de¬ 
generacies and interactions leads to interesting physics 
such as Mott phases with staggered loop currents [2^ 
1^ . Following the experiment of Atala et al. Q , Wei and 
Mueller [s^ used a variational approach to explore the 
effects of weak, short-range repulsive interactions, find¬ 
ing an additional phase, dubbed the biased ladder phase, 
where the density is uniform but different on each ladder, 
thereby breaking global Z 2 reflection symmetry. Here 
we generalize the theory and results of Ref. to long 
range dipolar interactions, finding a rich phase diagram 
as a function of the dipolar interaction strength, the syn¬ 
thetic magnetic field, and the relative tilt angle between 
the external field polarizing the dipoles, and the plane of 
the ladder. 

Our main results are summarized below: 

1. Generally long range dipolar interactions either de¬ 
stroy, or reduce the region of stability of the Meiss¬ 
ner or saturated chiral current phase. 

2. Repulsive dipolar interactions produce an inter leg 
charge density wave (GDW) phase where the densi¬ 
ties along the left and right legs of the ladder mod¬ 
ulate out of phase with one another. 

3. Next nearest neighbor interactions support a fully 
modulated biased ladder phase, where all the par¬ 
ticles are located on one leg of the ladder. 

4. Arbitrarily weak attractive interactions along the 
rungs destroy the Meissner phase entirely, and lead 
to a cascade of first order transitions between dis¬ 
tinct modulated density wave phases with different 
wave-vectors. 

5. Attractive nearest neighbor interactions along the 
ladder produce a regime of parameters where the 
biased ladder phase is the stable ground state at 
weak rung hopping. 


II. THE MODEL 

Our Hamiltonian for a two-dimensional spinless Bose 
gas on a two-leg ladder with lattice spacing a, takes the 
form (see Fig. [T]): 


n = no + nint (i) 

where Hq reads: 
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FIG. 1: (Color Online) Dipolar Bose gas on a two-leg 
ladder. The arrows on the sites of the ladder correspond to 
dipoles, aligned in the direction of the external field indicated 
by B. The external field could be an electric field in the case 
of polar molecules or a magnetic field for magnetic atoms such 
as Dy and Er. The hopping along the x direction is denoted 
by J. An artificial flux threads the system, which implies 
that hopping in the y direction picks up a position dependent 
phase The local on-site interaction is denoted by U and 
is always assumed to be repulsive U > 0. Additionally, there 
is a nearest neighbor interaction along x (ladder) and y (rung) 
denoted 14 and Vy respectively, and a next nearest neighbor 
interaction 14 nn. The sign and magnitude of 14, I 4 and 
Fnnn can be tuned by tilting the external field relative to the 
plane of the ladder. Higher order contributions to the dipolar 
interaction are significantly smaller in magnitude, and do not 
have much effect on the overall phase diagrams we present. 


where and ai^R are bosonic annihilation operators on 
the left (L) and right (R) legs of the ladder at position 
1. J and K denote the hopping matrix elements along 
the legs and rungs of the ladder respectively, and (j) is 
the magnetic flux per plaquette Q. This single particle 
model was introduced by Atala et al. Q. 

We further assume that there is an external (real) 
magnetic or electric field which polarizes the magnetic 
atoms/polar molecules, freezing out any internal degrees 
of freedom. As a result, the interaction Hamiltonian only 
includes density-density interactions and takes the form: 

^int = + (3) 

where fx, fi G (L^R), and U denotes the repulsive onsite 
interaction potential, which includes s-wave and dipolar 
contributions. The non-local part of the dipolar inter¬ 
action given by = V{1 - Scos^ )|r|-3, 

where V = D'^/a^ {D is the dipole moment), couples 
both left and right legs of the ladder. Here we define 
9i y ^ to be the angle between the external polarizing 
field and the vector r made by the lattice sites {I, /r} and 

Although the non-local dipolar interaction couples 
sites arbitrarily far apart, we restrict our calculations to 
nearest neighbor (NN) and next-nearest neighbor (NNN) 
interactions as shown in Fig.[TJ Previous studies on dipo¬ 
lar bosons in optical lattices have shown that retaining 
higher order terms has little qualitative effect on the over¬ 
all phase diagram [s^, but modifies the precise location 
of the phase boundaries. 
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The anisotropic nature of the dipolar interaction 
means that the interactions between sites on the same 
leg or the ^-direction) and sites 

on opposite legs aj or the y-direction), 

can in principle be different, and can be experimentally 
controlled by changing the alignment of the external field 
with respect to the plane of the ladder (see Fig. [T|). 
We therefore separately denote these interactions as 14 
and Vy, and denote the next-nearest neighbor interaction 
as yNNN- 

Generally, the on-site (U) and dipolar (V) interaction 
strengths are tunable by changing the transverse confine¬ 
ment The contact interaction is separately tunable 
using Feshbach resonances [H, [s^. We only consider 
repulsive on-site potentials in this work (^ . In experi¬ 
ments, the sign of the non-local part of the dipolar inter¬ 
action V is separately tunable from repulsive to partially 
attractive by tilting the external applied field at an ar¬ 
bitrary angle (0,x) relative to the ladder plane. For a 
general polar angle ff and azimuthal angle x, the interac¬ 
tion potentials along the ladder and rung direction read: 
14 PC V{l — 3 sin^ 9 cos^ x) a-nd Vy oc F(1 — 3 sin^ 9 sin^ x) 
1^ . A major advantage of this is that the physics of 
attractive ladders can be accessed, while ensuring over¬ 
all mechanical stability from a repulsive local interaction 


approach follows that of Wei and Mueller , who con¬ 
sidered the variational ground state wave-function for N 
particles, restricted to the lowest band: 

IGfeo) = + sin7/3lfeo)^|0) (4) 

where P±ko rire the annihilation operators for bosons at 
k = ±fco. The original boson operators can be expressed 
in terms of /3fc as akL = - sin ^/3fc and UkR = cos 

where the angle tan0fc = 2 sinksinci >/2 ■ Here |0) denotes 
the vacuum state, and 0 < 7 < 7 r /2 for feo > 0, and 7 = 0 
for k[j = 0. Absent interactions, this is the ground state 
for any 7 , but arbitrarily weak interactions will break 
this infinite degeneracy. 

B. Observables 

The local density on each leg is defined as = 
{Gko\a\^ai^y,\Gko)- We define the average density as 
n = A^/n, where G is the volume of the system. The 
expressions for the average de nsity on each leg and its 
modulations are derived in Ref. , but reproduced here 
for completeness: 


Completely attractive dipoles can also be obtained us¬ 
ing a time-dependent external field [4l|. In addition to 
the case where the external field is polarized perpendic¬ 
ular to the ladder, where all the interactions are repul¬ 
sive: Vx = Vy = V, Vnnn = Vl2\/2, we separately 
present the general phase diagram for attractive rung 
{Vy < 0) and ladder interactions (I 4 < 0) respectively. 
For convenience, we focus on two particular tilt config¬ 
urations: (i) external field along the rung (y) direction. 
Here Vy = -2V, 14 = V, Vnnn = and (ii) ex¬ 

ternal field along the ladder direction. Here I 4 = —2H, 
Vy = V and Vnnn = — H/4-\/2. Our qualitative results 
however are general, and robust to tilting the external 
field away from these angles, as long as the sign of the 
interaction along the ladder or rung direction does not 
change. 

The single-particle Hamiltonian in momentum space 
reads T-L{k) = —2Jcosfccos| -I- 2JsinA:sin|cr^ — 

Kux, and is readily diagonalized to yield two 
bands, whose energies are E± = — 2Jcosfccos(/)/2 ± 

4 sin^ k sin^ 4>l‘^ + [sB, IS IS ■ The lowest band 
has two minima a.t k = ±fco for K < 2 Jtand>/2 sind'/2, 
and a single minimum at fc = 0, for K greater than this 
value. 


A. Variational Approach I 

In this work, we consider two complimentary varia¬ 
tional approaches, valid for weak interactions, to study 
the ground state phase diagram of Eq. ([1]). The first 


ni,L = cos^ 7 sin^ 9kg + sin^ 7 cos^ 9ko + (5) 

- sin27sin0fep cos( 2 fco/) 
ni^R = sin^ 7 sin^ 9kg -f cos^ 7 cos^ 9kg + 

isin 27 sin 0 fcu cos( 2 A:oZ)- 

A key limitation of variational approach I is that the 
modulations on the L and R legs of the ladder are identi¬ 
cal. While this is a valid approximation for short ranged 
interactions considered in Ref. (S - where the interac¬ 
tion does not couple the two legs of the ladder, this is no 
longer the case for dipolar interactions. 

Working in momentum space, we define the current 

jfj. = {Gkg\J2k^l,y.^^k,y\Gkg). The total current jnet = 

jji is identically zero in equilibrium, but the chiral 
current jg = ~ jA can be finite. In order to visualize 

phases, it is also helpful to define the current in real space 
along the legs and rungs of the ladder as [^ : 

jL = - al^ai+i^y) ( 6 ) 

jt = iK{e~‘^'-‘^a] 

This is particularly useful for comparison with real space 
numerical methods such as DMRG. Explicit expressions 
for the rung and ladder currents in real space can be cal¬ 
culated using the variational wave-function Eq. and 
are provided in the Appendix. The net and chiral cur¬ 
rents are defined as j„et/c = J2iUIl ^ ^Ir) respectively. 

Using variational approach I for short range interac¬ 
tions {V = 0), Wei and Mueller [s^, found three phases: 
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1. Meissner (or saturated chiral current) phase. Here 
the rung current vanishes, and equal and opposite 
currents flow along the legs of the ladder. The total 
density is uniform. 

2. A biased ladder (BL) phase, where the density is 
uniform but different on the left and right legs, 
breaking global reflection symmetry, in addi¬ 
tion to the global C/(l) symmetry of the conden¬ 
sate. The rung current is zero. 

3. A modulated density (or vortex) wave phase, where 
the density is identical on the left and right legs but 
oscillates with a period incommensurate with the 
underlying lattice. This phase breaks additional 
U{1) symmetry, associated with translations of the 
density wave. 

The presence of local currents on each plaquette in the 
modulated density phase are analogous to vortices in a 
type-II superconductor. The terminology of vortex and 
Meissner ^ases was first introduced by Orignac and Gi- 
amarchi [2^ . and Atala et al. Q interpreted their data 
using this language. Here we will interchangeably use 
the terminology saturated chiral current (or Meissner) 
and modulated density wave (or vortex) phase. We cau¬ 
tion however that for strong interactions, vortex states 
may not always be accompanied by density modulations, 
indeed homogeneous insulating vortex phases have been 


found in DMRG studies [s^. As we explicitly derive in 
the Appendix, in all the three phases, the net ladder cur¬ 
rent jnet remains zero on every plaquette. 

We note that in principle, one can construct general 
non-local density and bond operators such as 

for sites i,j and £ {L,R} [i2,li3i which can take 
on nonzero values, especially for systems with long range 
interactions. In the weakly interacting limit however, 
whenever {Pko) becomes finite, these other order param¬ 
eters automatically take on non-zero expectation values. 
Therefore it suflices to only consider the density here. 
How these non-local order parameters vanish near the 
Mott transition however, is a separate question. 

The ansatz of Eq. o can be readily employed to cal¬ 
culate the total energy of the dipolar gas E = E-{k) + 
{Gko\'Hmt\Gko) 1 where we treat 7 and fco as variational 
parameters. The saturated chiral current (biased ladder) 
phases correspond to 7 = 0 and fc = 0 (y^ 0 ) respec¬ 
tively, whereas the vortex or modulated density phase 
corresponds to 7 0, fc 0. A new feature of dipolar 

systems is that 7 can be pinned to a finite value, even 
as fc —>■ 0 , which implies a state with long wave-length 
density modulations, yet a fully saturated chiral current. 
This state breaks the same symmetries as the modulated 
density-wave phase. 

The variational form for the interaction energy is: 


A(fco,7) = E-{ko) + 


[U + 14)n 


1 


1 


1 


(1 - - sin"' 9ko - 2 + 2 27 sin"' 9ko) + 


([/-I -14 cos(2A:o))n , 2 


(14 -I- 2Hnnn)?^ , . 2n ,■20 ■ 2n •2o\i (^y + 2Hnnn COS (2A:o))n ^._^2 


l(sin'^ 9ko + sin'^ 27 - sin'^ 9kg sin^ 27) -|- 

I- 


16 


sin'"' 27 sin'^ 
sin'"' 27 sin^ 9ko 


One readily checks that in the absence of non-local dipo¬ 
lar interactions, Vx = Vy = Vnnn = 0 , the interaction 
energy reduces to the form obtained in Ref. [s^. 


C. Variational approach II 

It is useful to note that the three phases found in 
Ref. [15 are directly analogous to the corresponding 
phases of a weakly interacting spin-orbit coupled Bose 
gas in the continuum [43| . which has been extensi vely 
studied theoretically and experimentally [ 2 .l4.flOl ldSl - H^ . 
There, the one-dimensional spin-orbit coupling plays the 
role of the artificial magnetic field. The modulated den¬ 
sity wave phase in the ladder problem corresponds to a 
stripe phase of the spin-orbit coupled Bose gas, where 
both single particle minima are occupied; the biased lad¬ 
der corresponds to the situation where only one mini¬ 
mum is occupied, and the Meissner (saturated chiral cur¬ 


rent) phase corresponds to a non-magnetic phase, where 
ko = 0 . 

Motivated by this, we consider a second variational 
ansatz, first introduced by Li et al. [15 to study the spin- 
orbit coupled gas, but adapted to the present problem. 
This has the key advantage that it allows us to readily 
generalize to a situation where the modulations on the 
left and right legs of the ladder are unequal. We replace 
the boson operators on the left and right legs of the ladder 
by classical fields: 

ai,L = sin + G 2 L cos ( 8 ) 

ai,R = 

where Ciy,,C 2 ^ are complex numbers which satisfy the 
constraint: sin^ ^{\Gil\^ + \G 2 r\^) + cos^ %^(|C' 2 lP + 
IGifip) = 1, from number conservation. Note that the 
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variational approach of Refs. [s^ . 1^ corresponds to the Inserting Eq. ([81) into the Hamiltonian, the expression 
limit where Cil = Cm and C 2 L = C 2 r- This restricted for the total energy reads: 
ansatz yields a density wave texture, but the spin den¬ 
sity remains homogeneous. However, in the presence of 
inter-rung interactions, our more general ansatz is natu¬ 
rally able to capture states with spin and density wave 
modulations, which are given by n; = ni^L + ni,R and 
Si = nim - ni,L where: 


n^L = ICiLpsin^ ^ -b \C2l? cos^ ^ + (9) 

^(C'iLC2\e2*'=o'sin4o+c.c) 
ni,R = |C2flP sin^ ^ -b \Cir\^ cos^ ^ -b 
i(Ci«C2*fle2*'=“'sin0fe„+c.c) 


Emf = E_{ko) + (dCiil" + \C2r\^) sin^ ^ + \C2l\^) cos^ + (10) 


^feo 


f^ko 


(U + 14 COS^ fco)n 


sin2 0fco(|CiLnC2Lp + |CiijnC2flP) ■ 


V,,n 


{\Cil\'^\C2r\ 


^sin4^ + |C2LnCiflpcos4^ 


,4 

2 


i sin^ 9 k„i\CiL\^\CiR\^ + \C2 l\^\C2r\^ + C*mC*2RC2LCiR + c.c)) + VNNNn{\CiR\^\C2L\^ cos 

\Cil\'^\C 2 r\^ sin^ ^ + l sin2 ek„i\CiL\^\CiR\^ + \C 2 l\^\C 2 r\^ + C*mC* 2 RC 2 LCiRCOs 2 ko + c.c)) 


Note that unlike the total energy of variational ap¬ 
proach I, Emf is no longer symmetric under k^ —> — fco- 
This is because under k^ —>• — feg, cos —>■ sin^. 

Therefore, for generic values of the co-efficients and 
C 2 /i, EMF(fco) ^ Emf(— fco)- One way to impose this 
symmetry is by setting Cil = Cir and C 2 L = C 2 R in 
the variational ansatz. Here we look for more general so¬ 
lutions by numerically minimizing the variational energy 
Emf for arbitrary complex and € 2 ^., with the con¬ 
straints of particle number conservation and EMF(fco) = 
Emf(— fco)- 

In what follows, we use both variational approaches to 
obtain the ground state energy. For repulsive interac¬ 
tions along the rungs, the inter leg charge density wave 
solution found using variational approach H always has 
lower energy than the modulated density wave phase of 
variational approach I. All other phases we find are cap¬ 
tured by both ansatzes. Moreover, both approaches yield 
identical phase boundaries. 

Before turning to the results, we comment on the va¬ 
lidity of variational, mean-field approaches in ID, where 
fluctuations are important and destroy long range or¬ 
der. The crucial parameter controlling the validity of 
our approximations is the ratio of the interaction to the 
kinetic energy C = Eint/Ekin. In ID, the interaction en¬ 


ergy scales as Eint Un, while the kinetic energy scales 
as Ekin ^ /2mcP, where d is the mean interparticle¬ 
spacing d oc 1/n. The weakly interacting, mean-field 
regime occurs when the parameter C oc E/n <C 1 [i^ . 
Thus our approximation works best at high densities, 
which corresponds to having a large number of bosons 
per site, which is precisely the regime we consider here. 


III. PURELY REPULSIVE DIPOLES 
A. Qualitative Features 

We begin by examining the phase diagram of the two- 
leg ladder for purely repulsive dipoles. Before turning to 
the numerical results, we discuss the qualitative physics 
we expect from long range interactions. We first consider 
the role of the ladder and rung interactions separately, 
and then present the full phase diagram. 

The nearest neighbor interaction [V^ni^^ni+i^) along 
the ladder can be eliminated by placing particles on either 
the even or odd sites of the lattice, which corresponds to 
a modulated density wave phase. By contrast, it penal¬ 
izes states with homogeneous density profiles on the legs 
of the ladder, which is the Meissner and biased ladder 
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FIG. 2: (Color Online) Repulsive Dipoles: (a) Global phase diagram showing the density difference between left and right 
legs of the ladder as a function oiVj J and Kj J aXU = Q and (j) = -k 12. Dipolar interactions push the Meissner phase to stronger 
rung-ladder coupling and give rise to an interleg CDW phase, where the relative densities on each leg of the ladder modulate 
out of phase with one another in real space as shown in (c). At intermediate rung hopping and weak dipolar interactions, a 
biased ladder (BL) phase is present, (b) Real space density profiles in the left (solid), right (dashed) and total density (dotted) 
for large next nearest neighbor interactions. All the particles reside on only one leg of the ladder, producing a fully modulated 
biased ladder phase, (d) Real space current profile at <() = O.Ott, showing switching of the sign of the chiral current. Dhar et al. 
refer to this phase as a chiral superfluid, (e) Chiral current at ()> = O.Qtt, [7 = 0, V/J = 1 and K/J — 1.5. 


phases. We thus expect the dipolar ladder interaction to 
suppress these phases at large rung to ladder coupling 
strength (K/J). 

To see this, note that in Eq. o. this interaction con¬ 
tributes to the on-site interaction, and also yields a mo¬ 
mentum dependent term term proportional to cos(2fco), 
which can compete with the local interaction U. In the 
following, we set Ennn and Vy to zero in Eq. m for sim¬ 
plicity. Although this choice is somewhat artificial, it 
is instructive in developing a systematic understanding 
about the effect of long range interactions. 

For purely contact interactions, when K > K^, the 
single particle minimum occurs at fc = 0, and the en¬ 
ergy of the Meissner state (fc = 7 = 0) reads: E_(0) -I- 
17/4 < if_(fco) -I- 17/4 -I- 17/8 sin^ 0fc(,, which is the en¬ 
ergy of the modulated density state at fc = fco and 
7 = 7r/4. Thus the Meissner phase wins at large K. 
By contrast, for nearest neighbor dipolar interactions 
along X, the energy of the modulated density state reads: 
i?_(fco)-l-t4/4-|-I4 cos(2fco) sin^ which can be made 
lower than the Meissner phase energy (if_(0) -|- 14/4), 
by choosing 7r/4 < fco < tt/2, provided 14 is sufficiently 
large. Therefore, for any value of K, however large, there 
is always a transition out of the Meissner phase into a 
modulated density wave phase, provided I 4 is made suf¬ 
ficiently large. For typical parameters K/J ^ 1.5, the 
Meissner phase is fully destroyed in favor of a modulated 
density wave phase for 14/J 7. 

We now consider the effect of a repulsive inter-rung 
interaction on the phase diagram, which will naturally 
be present, whenever the external field is polarized per¬ 
pendicular to the plane of the ladder. While this term 
also favors non-zero fco and 7 = 7r/4, namely a modu¬ 
lated density phase, repulsive rung interactions penalize 
the left and right rungs from having identical density 


modulations. Indeed, we would expect that strong inter¬ 
rung repulsion would lead to “phase separation”, where 
the density modulations on the left and right legs are 
out-of-phase with one another. We term such a state an 
interleg charge density wave (CDW). This state therefore 
breaks local Z 2 reflection symmetry, even though glob¬ 
ally, Z 2 reflection symmetry is unbroken, as the densities 
oscillate about the same average value. In the absence 
of gauge fields, such a state would be the ladder analog 
of a checkerboard solid phase of dipolar bosons in optical 
lattices [13, However, our state additionally has non¬ 
zero rung and ladder currents arising from the synthetic 
flux threading the ladder. We remark that the interleg 
CDW phase we find here is the weak coupling an alog of 
the “SDW” phase found by Petrescu and Le Hur [3l| in 
the strong coupling limit of our model (absent next near¬ 
est neighbor interactions) at unit filling fraction. Fur¬ 
thermore, unlike the checkerboard phase of dipolar lat¬ 
tice bosons or the SDW phase, the interleg CDW phase 
we find is generally incommensurate with the underlying 
lattice, as the wave-vector of the density wave oscillations 
is set by the minima in the single particle dispersion. 

B. Inter leg CDW 

The global phase diagram found by numerically mini¬ 
mizing the total energy Fmf (Eq. 10) for I 4 = I 4 = H 
and U = Vnnn = 0, and (j) = 7r/2, is shown in Fig. [5Ka). 
Indeed for repulsive dipolar interactions, we find that a 
large portion of the phase diagram is occupied by the in¬ 
terleg charge density wave phase. From Eqs. we plot 
the densities on the left (solid) and right (dashed) legs 
of the ladder, which are oscillating out of phase with one 
another. 
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As expected from our discussion above, the Meissner 
phase indeed gets pushed to larger Kj J with increasing 
dipolar interactions. We have checked that variational 
approach I yields the same qualitative phase diagram and 
identical phase boundaries, but that theory cannot dis¬ 
tinguish between a modulated density wave phase and 
an inter leg CDW phase. The latter always has lower en¬ 
ergy for repulsive rung interactions. For even larger V 
(not shown), the Meissner phase disappears completely, 
as we argued from our qualitative discussion above. The 
transition from the interleg CDW to the Meissner phase 
is always first order for repulsive dipoles. For interme¬ 
diate K, and weak dipolar interactions there is a biased 
ladder phase, which spontaneously breaks global Z 2 re¬ 
flection symmetry, as the average densities on the left 
and right legs are no longer equal. For very weak dipolar 
interactions, our numerics find a modulated biased lad¬ 
der phase, where the densities on each leg modulate in 
real space but about different average values. However, 
if such a phase exists, we believe that it should only be 
present in a small window of the phase diagram. 

Inclusion of the next nearest neighbor term I^jnn, 
has no qualitative effect on the phase diagram presented 
above, as the next nearest neighbor interactions are sig¬ 
nificantly weaker than the ladder and rung interactions. 
Quantitatively however, next nearest neighbor interac¬ 
tions increase the window of stability of the biased ladder 
phase. Interestingly, if the next nearest neighbor inter¬ 
action is the largest interaction in the problem, then the 
ground state is a fully modulated biased ladder, where 
all the bosons reside on only one leg of the ladder (see 
Fig. [2Kb)). This state not only breaks global Z 2 re¬ 
flection symmetry, but also has an additional Goldstone 
mode associated with breaking translational symmetry. 
This state is the ladder analog of the stripe phase of 
two-dimensional lattice dipolar bosons. Although such a 
phase is theoretically interesting, this parameter regime 
cannot be attained with dipolar interactions alone. 

Tilting the dipoles away from the plane perpendicular 
to the ladder does not affect any of the qualitative con¬ 
clusions of Fig. 121 provided the nearest and next nearest 
neighbor interactions remain repulsive. The precise loca¬ 
tions of the phase boundaries however will shift. 

Dipolar systems have long been sought after as ideal 
candidate systems for exploring “supersolidity”: a phase 
of matter which exhibits dissipationless flow analogous to 
a superfluid while simultaneously possessing crystalline 
order (50j |. However a major experimental challenge in 
observing supersolidity in dipolar systems is that the 
dipolar interaction strengths in most magnetic atoms and 
polar molecules that have been cooled to date are far too 
weak, and are easily overwhelmed by the contact inter¬ 
action, which produces a homogeneous superfluid [s^ . 
However, in the present system, translational symmetry 
breaking is not an interaction effect, rather a single par¬ 
ticle one, owing to the multiple minima in the single par¬ 
ticle dispersion. Having broken translational symmetry 
to produce a modulated density wave phase, there is no 


additional on-site energy cost to displacing the density 
modulations on the left and right legs of the ladder. But 
this lowers the overall dipolar energy by avoiding the rung 
interaction. Therefore a checkerboard state or an interleg 
CDW phase can occur for very weak dipolar strengths in 
this system, and survives even in the presence of a non¬ 
zero U. 


C. Chiral current switching 

We briefly discuss what happens at other values of 
the flux. Although the experiments do not tune the 
flux directly, DMRG studies have extensively explored 
the phase diagram for different values of the magnetic 
field [H, [ 13 . Remarkably, we find that near (but not 
equal to) </> = tt, the sign of the chiral current develops a 
switching pattern from plaquette to plaquette as shown 
in Fig.|2K(d) and (e)). In Fig.|2Ke), we plot the chiral cur¬ 
rent ai (j) = O.Qtt, which shows this sign reversal. A sim¬ 
ilar pattern is predicted to occur in the fully frustrated 
Bose Hubbard model studied by Dhar et al. [^, , who 

refer to this state as a chiral superfluid. We remark that 
this is different from the findings of Greschner et al. [s^ , 
where the global chiral current reverses sign compared 
to the chiral current at small fluxes. Here we find that 
the sign of the global chiral current is still the same as 
that for small fluxes. We find that this switching pattern 
occurs even for weak interactions, at large enough values 
of (p, when the Meissner phase is destroyed in favor of 
a vortex phase. Increasing K increases the magnitude 
of the chiral current at fixed flux, and should make this 
effect experimentally observable. 

IV. PARTIALLY ATTRACTIVE DIPOLES 

We now turn to the study of partially attractive 
dipoles, which can be realized by tilting the external field 
relative to the plane of the ladder. For concreteness, we 
focus on two particular cases: (i) external field aligned 
parallel to the rungs (attractive rung and next nearest 
neighbor interaction) and (ii) external field aligned par¬ 
allel to the legs of the ladder (attractive ladder and next 
nearest neighbor interactions). We emphasize that even 
though we focus on specific tilt angles below, our results 
in Sec IV A and B are valid over a wide range of tilt 
angles as long as I 4 <0,14 >0 and I 4 < 0, Hy > 0 
respectively. As before we calculate the phase diagram 
using both variational approaches. Where comparison is 
possible, they yield identical results. 

A. Attractive Rung interaction 

We first consider the case where the rung interaction 
Vy is attractive and equals Vy = —2V, while the interac¬ 
tion along the ladder remains repulsive Vx = V. We also 
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FIG. 3: (Color Online) Dipoles pointing along y (rung) direction: (a): Phase diagram of the attractive rung, dipolar 
bosonic ladder for [7 = 0, and 0 = 7 r /2 as a function of V and K/J. A modulated density wave phase (with 7 = ^r/d) is 
found throughout the entire phase diagram. Density plot shows evolution of fco with V/J and K/J, which approaches zero 
as K becomes large, but 7 remains finite. Inset shows real space total density at two different values of K, at fixed V. ((b) 
and (c): Evolution of the chiral current and fco respectively for fixed V, as a function of K/J, showing jumps. The discrete 
jumps indicate first order transitions between CDW phases with different incommensurate wavelengths and is reminiscent of 
the devil’s staircase pattern found in dipolar Mott insulators and superfluids in ID [T^, . For comparison the evolution of fc 

for purely short range interactions is also shown (dashed), revealing a single first order transition from the modulated density 
phase to the Meissner phase, (d) Schematic showing the real space density and current profile in the two-leg ladder. Increasing 
K/J increases the wavelength of the modulations (shaded region). Arrows indicate the direction of the chiral current In this 
picture, “vortices” sit in the regions between the shaded areas, where the chiral current reverses direction. Therefore, the 
vortex density decreases as the K/J increases (c.f. Ref. [33] Fig. 5). Both variational approaches never find a transition to a 
homogeneous density (Meissner) phase at any K. 


turn on a finite, attractive NNN interaction. As in the 
repulsive case, we consider a purely dipolar interaction 
(7 = 0. Owing to the attractive rung interactions, the 
density modulations on both legs are in phase with one 
another, and both variational approaches yield a modu¬ 
lated density wave phase as the ground state. 

In Fig. EKa), we plot the phase diagram of the two 
leg ladder, as a function of K/J and V/J obtained from 
Eq. ([7]). The density plot reveals that a vortex or mod¬ 
ulated density phase (7 = 7 r/ 4 ) occurs throughout the 
entire phase diagram, but that the wave-length of the 
modulations (or the vortex core radius) diverges as K/J 
is increased. At large K/J, the chiral current saturates 
to its maximal value (fco 0 ), but 7 remains pinned to 
7 r/ 4 , so the density modulates in space. Thus we con¬ 
clude that arbitrarily weak attractive rung interactions 
destroy the Meissner phase. 

To understand the destruction of the Meissner phase, 
we once again consider the energy difference between the 
Meissner and a modulated density wave phase. Set¬ 
ting Vnnn = 0, which has no qualitative impact on 
the physics, the energy of the Meissner phase from 
Eq. ([7]) reads: E = £’_(0) — E/8 while that of 
the modulated density wave phase is E-{ko) — E/8 -I- 
Esin^(dfc(,)/8(l/2cos(2fco) — 1). At large K, where the 
Meissner phase usually occurs, the energy difference be¬ 
tween E- {0)—E- (fco) approaches a constant independent 
of K, E-{0) — E-{ko) —>• —•\/2(l — cos(fco)) . Therefore, 
for a modulated density wave phase to win, we require: 
E sir?{9ko)/^{1/2 cos(2fco) — 1) < —•\/2(l — cos(fco)). It 


is easy to see that this condition can be satisfied even for 
arbitrarily weak E, as long as fco —f 0. Therefore, for any 
value of E, there is a critical fccrit, such that for fco < fccrit, 
the modulated density wave phase with wave-vector fco 
and 7 = 7 r /4 has lower energy than the Meissner phase. 
Hence the ground state is always a modulated density 
wave or vortex phase. 

Turning on a finite U however, restores the Meissner 
phase. Indeed, it is easy to show from Eq. that for 
fco —>■ 0, a hnite 7 occurs only when U < Uc = V/J\/2. 
This condition can achieved experimentally by using Fes- 
hbach resonances to tune U near a zero crossing. We 
expect that including higher order terms will increase 
Uc, as longer ranged interactions will be more attractive. 
However, we expect the quantitative corrections to be 
small, as longer range interactions decrease in magnitude 
as 1/r^. 

Physically the destruction of the Meissner phase for 
attractive rung interactions can be simply understood: 
Absent short-range repulsive forces {U = 0), the dipolar 
interaction favors the formation of a CDW, and the at¬ 
tractive rung interaction implies that the CDWs on both 
legs will be in-phase with one another, i.e., a modulated 
density wave phase. Upon increasing K, the two single 
particle minima at ±fco move closer to one another, but 
the total dipolar interaction energy can be lowered by 
condensing at both minima, rather than only in a single 
minimum, which is the case for repulsive short-range in¬ 
teractions, which favors a homogeneous density. When U 
becomes larger than a critical strength, the homogeneous 
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density Meissner phase is restored. 

Strikingly, the density plot reveals that feg shows a se¬ 
quence of jumps as it varies from = 7r/4 to 0 with 
increasing K (Fig. EJc)). These jumps are also mani¬ 
fested in the chiral current (Fig. [3Kb)), and are indicative 
of a cascade of first order transitions between modulated 
density wave phases with different wave-vectors. The 
series of first order transitions between different mod¬ 
ulated density wave phases are reminiscent of the Lut- 
tinger staircase feature observed in quasi-lD dipolar su¬ 
perfluids [l^. There, the dipolar Luttinger liquid be¬ 
comes unstable towards the formation of a cascade of 
solids with fractional filling, in the presence of an arbi¬ 
trarily weak lattice potential. Although we work in the 
opposite limit (namely of large incommensurate occupa¬ 
tion), our theory nonetheless yields a similar cascade of 
CDW phases, where the periodicity of the density wave is 
tuned not by filling fraction, but rather by the hopping 
parameter. For comparison, in Fig. EKc) we also show 
the evolution of fco for a gas with contact interactions U 
of the same magnitude, which shows a single first order 
transition from a modulated density wave to a Meissner 
phase. 


B. Attractive Ladder interaction 
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FIG. 4: (Color Online) Dipoles pointing along x (lad¬ 
der) direction: Density plot showing the density difference 
between the two legs of the ladder normalized to the total den¬ 
sity for attractive ladder interactions: lA = —2V, Vy = V, 
and U/J = 1. For finite V/J, there is a transition from an 
interleg CDW to a biased ladder phase at arbitrarily small 
K/J, which occupies a large portion of the phase diagram. 
Increasing K/J leads to a direct transition from the biased 
ladder to the Meissner phase. Increasing the dipolar interac¬ 
tions further leads to a collapse of the gas, signaled by the 
appearance of imaginary frequencies at long wave-lengths. 


We now turn to the case where the magnetic field is 
tilted along the ladder (x) direction. We therefore have 
Vx = —2F, Vy = V. Numerically minimizing the varia¬ 
tional energy Eq. (10), the ground phase diagram for this 
case is shown in Fig. |4| We have set U/J = 1. 

For weak ladder interactions V « U, the phase dia¬ 
gram resembles that of a gas with local interactions: a 
direct first order transition from an interleg CDW phase 
to a Meissner phase is found. (Note that for the choice of 
U, there is no intervening biased ladder phase (Ref. [s^l.) 
However, increasing attractive dipolar interactions de¬ 
stroys the interleg CDW phase in favor of a biased lad¬ 
der phase for arbitrarily weak rung to ladder couplings. 
Increasing K/J leads to a second order transition from 
the biased ladder phase to the Meissner phase. We re¬ 
mark that the biased ladder phase has not yet been ob¬ 
served ex per imentally, despite being predicted theoreti¬ 
cally [s^, 13^ . [blj . The tilted dipolar system displays a 
wide biased ladder regime, and may be ideally suited for 
observing this phase experimentally. 

A similar phase diagram, where the interleg CDW 
phase is replaced by a modulated density wave phase 
is expected for purely attractive short range interactions 
U < 0. However, such a gas is not mechanically stable 
at zero t emp erature, and collapses at sufficiently high 
densities [^,[^. In our case, the appearance of a biased 
ladder phase at small values of the rung hopping is driven 
entirely by attractive long range forces, and the gas can 
still be stabilized by maintaining a repulsive on-site po¬ 
tential. 

To compute the stability of the biased ladder phase. 


we calculate the excitation spectrum about the biased 
ladder ground state \Gko) = l/-\/M(/3^^)^|0). Writing 
fluctuations above the ground state as /3k = VNSk,ko + 
(1 ~ 5k,ko)lk-ko [13 ! '^6 expand the Hamiltonian Eq. ([T]) 
to quadratic order in the 7 operators. Diagonalizing 
the effective Bogoliubov Hamiltonian H — fiN, where 
^ = E_{ko) + {U + Vx)n{sm'^ (^'fco/2) + cos'^ (Oko/^)) + 
{Vy + VNNN/’\/ 2 )?^sin^ ( 0 fco/ 2 ) cos^ ( 0 fe(,/ 2 ) is the chemi¬ 
cal potential, yields the low energy excitation spectrum. 
As in the case of short-range interactions, the spectrum 
has a characteristic roton-maxon feature [36j |. which has 
also been recently found in the shaken lattice experiment 
of Parker et al. The dipolar interaction pushes 

the roton minimum to larger k. Near the interleg CDW 
to biased ladder phase transition, the roton minimum is 
present for arbitrarily small K, and disappears as K is 
increased. 

The instability towards collapse is signaled by the ap¬ 
pearance of imaginary frequencies at long wave-lengths. 
For U/J = 1, the gas becomes dynamically unstable to¬ 
wards collapse at V/J ~ 0.5, or I 4 ^ —U. The critical 
interaction strength at which collapse occurs is weakly 
dependent on K. Increasing K/J pushes the collapse to 
somewhat stronger interactions. 

V. SUMMARY AND DISCUSSION 

Experimentally, our system can be realized using 
highly ma gne tic dipolar atoms, Rydberg atoms or polar 
molecules (lll - ll^ . Magnetic atoms are ideally suited for 
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this study, as we do not incorporate spin degrees of free¬ 
dom, and therefore do not have any dipolar loss. More¬ 
over, coupling magnetic atoms to Raman lasers is be¬ 
lieved to have much lower heating losses [s^, which is a 
major advantage over conventional alkalis. With polar 
molecules, a key experimental challenge that needs to be 
overcome is the low densities needed to prevent chemi¬ 
cal reactions. Non-reactive molecules or quantum-Zeno 
methods of loss suppression , may be required in or¬ 
der to boost the densities to values suitable for studying 
many-body physics. 

All the phases discussed here can be probed experi¬ 
mentally by measuring the in situ density and the lo¬ 
cal currents Q. The interleg charge density wave phase 
breaks local Z 2 symmetry spontaneously. However, as 
the dipolar interaction is long ranged, even if there is no 
hopping between neighboring ladders, long range inter¬ 
actions will couple neighboring ladders, giving rise to an 
overall checkerboard pattern, which can be detected us¬ 
ing Bragg spectroscopy [53|. The modulated biased lad¬ 
der phase is more challenging to detect, as the density 
asymmetry between the left and right legs of the ladder 
will randomly change from ladder to ladder, averaging to 
zero. However, by detuning the Raman lasers from res¬ 
onance, an external bias can be applied [s^, which then 
polarizes all the ladders in the same way. 

We comment briefly on the relationship between our 
work and recent experiments on bosons and fermions in 
synthetic dimensions [1, . The physical rung in the lad¬ 

der we consider, corresponds to the synthetic dimension, 
implied by the internal spin states of the atoms. There¬ 
fore the onsite density-density and spin-spin interactions 
can be interpreted as effective “long range” interactions 
in the synthetic (rung) dimension. A crucial difference 
however is that in these experiments, long range interac¬ 
tions only occur in the synthetic direction; in the physical 
direction, there is a one-dimensional optical lattice and 
the interactions are still short-ranged. If the analog of 
an interleg CDW occurs in these experiments, it will be 
manifest as a spin-density wave, where the local spin ori¬ 
entation changes from site to site, such that on a given 
rung of the ladder, only one spin component is present. 

To conclude, we have used two complementary mean- 
field approaches to explore the interplay between large 
magnetic fields and long range density-density interac¬ 
tions on a two-leg bosonic ladder, finding a rich phase 
diagram. Our approximations are justified in the large 
density limit of large numbers of bosons per site, which 
for ID systems corresponds to the mean-field regime. As 
in the experiment of Atala et al. Q, we largely fix the 
magnetic flux, but modulate the rung-to-ladder hopping 
and the interactions. Generally we have shown that dipo¬ 
lar interactions destroy the Meissner phase completely or 
reduce its regime of stability. We expect that more so¬ 
phisticated bosonization treatments should yield similar 
qualitative answers, as it is already well known that dipo¬ 
lar Luttinger liquids are unstable towards CDW forma¬ 
tion [ 1 ^. A full bosonization treatment of this problem 


remains an exciting topic for further study. Repulsive 
density-density dipolar interactions lead to an interleg 
CDW phase, where the total density is uniform, but the 
densities along the left and right legs of the ladder modu¬ 
late in space, out of phase with one another. This phase 
is stable for weak next nearest neighbor interactions, but 
for strong next nearest neighbor repulsion, a fully mod¬ 
ulated biased ladder phase is found, where all the atoms 
are either on the left or the right leg of the ladder. For 
values of the flux near ^ tt, we obtain a modulated 
density wave phase where the chiral current switches sign 
from plaquette to plaquette. 

For purely attractive rung interactions, we find that 
the Meissner phase is completely destroyed. Instead, we 
find a cascade of first order transitions between CDW 
phases with different wave-vectors, reminiscent of the 
Luttinger staircase M- Whether this pattern develops 
into a full Luttinger staircase at low filling fractions will 
be studied in future work. Importantly, as these jumps 
can be tuned by changing the Raman coupling, this may 
open the possibility of observing this staircase pattern in 
experiments. 

Finally, we emphasize that although we have only con¬ 
sidered specific tilt angles in this work for simplicity, our 
results are more general in that the phase diagrams ob¬ 
tained here will survive small deviations away from the 
tilt angles we consider. For example, the biased lad¬ 
der phase is the ground state as long as 14 < 0 and 
1 14 1 < U and 14 > 0. This condition is met for all angles 
sin“^ I/a/S < 6 < TT12. Similarly, we expect an inter leg 
CDW ground state for small K/ J whenever Vy > 0. 
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Appendix A: Expressions for local ladder and rung 
currents 

To calculate the currents in Eq. we first 

transform our variables into Fourier space as fol¬ 
lows akL = l/VTiX); and auR = 

l/-\/DX]i Inserting these expressions into 
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Eq. (( 6 ]) and expressing in terms of /3fc, we obtain: 

K 

= — — sin27sin2fco^cos0fc(, (Al) 

The ladder currents are similarly found to be: 


A = ^ (2 sin^ ^ sin (fco + ^) cos^ 7 + 

2 cos^ sin (—/cq + sin^ 7 + 
j SingfcoSin7 ^^_2.fco;^^-.(fco+^) _ _ 


(A2) 


The reader may readily check that using the fact that 
there is a symmetry from fco —> — fcoi a-nd that cos^ = 
sin^ the net ladder current vanishes on every 

site I, and the net chiral current is jc = Thii.j'iL ~ Jib) ~ 
4Jsin^ ^ sin (feg + (/)/ 2 ). 


Jib = ^(-2 cos^ ^ sin (-fco + cos^ 7 - 

2 sin^ sin (fco + sin^ 7 + 
^ sin^fcp sin7 ^_2ikoi^^-i{ko-i) _ g-i(feo+|)^ _ 
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